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1 Introduction 



The idea of the dynamical symmetry breaking in quantum field theory has been in- 
troduced quite long ago JTJ. In the attempts to realize this idea in the early universe 
one should study the external gravitational field (or curved spacetime). The investiga- 
tion of the dynamical symmetry breaking in four-fermion models [p], 0] in four and three 
dimensional curved spacetime has been started in Refs. |3], f|] using the weak curvature ap- 
proximation (for a review and list of references, see ||). A validity of the weak curvature 
approximation is discussed in Ref. 0. In the case when the external gravitational field 
is treated exactly (de Sitter or anti-de Sitter or Einstein universe background) the phase 
structure of the four-fermion models has been studied in Refs.JTJ. (for a renormalization 
group approach, see also 0). 

There appeared recently some indications that early universe may contain large pri- 
mordial magnetic fields. The role of these fields in cosmology (in particularly, in the 
inflationary models of the universe) has been discussed in Refs.[f|. The presence of 
magnetic field in the early universe may lead to different effects. In particular, it is rea- 
sonable then to investigate the dynamical symmetry breaking in curved spacetime with 
magnetic fields^. Such study for 3D and 4D four-fermion models has been undertaken 
in Refs. [jl3| , |14|1 in the approximation when the effects of the gravitational and magnetic 
field maybe simply summarized (in other words, the correspondent terms in the effective 
potential are evaluated with the neglecting of the interaction between gravitational and 
magnetic fields). However, the coherent effect of combined gravitational and magnetic 
field may become very relevant. 

That is the purpose of the present work to investigate the phase structure of four- 
fermion models in weakly curved spacetime with constant magnetic field. We develop 
the approximation where the effective potential and Green function (GF) of the theory 
maybe represented as the expansion in the powers of curvature invariants. In each order 
of this expansion the external magnetic field is treated exactly, i.e. we take into account 
the coherent effect of combined gravitational and magnetic fields. 

The paper is organized as follows. In the next section we present new form of the local 
momentum representation of propagators [[Tj|] with account of the constant magnetic 
field. In each order on the expansion in curvature invariants for Green functions the 
magnetic field is included exactly. The detailed calculation is presented for scalar and 
spinor Green functions. Section 3 is devoted to the study of phase structure in the NJL 
model in curved spacetime with magnetic field. We calculate the effective potential with 
account of linear curvature terms for 3D and 4D models. The numerical investigation 
of the effective potential for 4D NJL model is given. In section 4 we discuss the SUSY 
NJL model non-minimally interacting with the external gravitational field and minimally 
interacting with constant magnetic field. The effective potential is evaluated and the 
correspondent phase structure is given. Some discussion and outlook are presented in 



§For the works where four-fermion models have been discussed in the presence of constant electro- 



magnetic fields, see lid, 11, [12 
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the last section. 



2 Scalar and spinor Green functions in the magnetic 
field in curved spacetime 

One of the fundamental objects in the field theory is a Green function. It is a basic 
quantity for a number of calculations of quantum effects. Here we calculate the scalar 
and spinor Green functions in the magnetic field in curved spacetime to study the phase 
structure of the NJL and SUSY NJL model in the external gravitational and constant 
magnetic fields. 



2.1 Green function for a scalar field 

First we study GF for the scalar field G(x,y;m) in an external electromagnetic and 
gravitational field. We assume that the spacetime curves slowly and neglect the terms 
involving the metric derivatives higher than third (weak curvature approximation). How- 
ever we treat the external magnetic field exactly. The GF for the scalar field G(x, y; m) 
satisfies the Klein-Gordon equation : 

(iDHDp - m 2 - £R)G(x, y; m) = -^8 D {x - y), (1) 



where £ represents the non-minimal coupling constant with external gravitational field(see 
p4|), m is the scalar field mass and D^ 1 = V M — ieA^. In the present paper the vector 
potential of the external electromagnetic field is chosen in the form 

M*) = -\f^\ (2) 

where is constant matrix of electromagnetic field strength tensor. The covariant 
derivative for a scalar field is only an ordinary derivative V^ 1 = d^. For a vector field t> M 
the covariant derivative is rewritten as 

v,/ = c^ + rv^ 

' dpy/=gv*. (3) 



Thus the first term of the left hand side in Eq.(|I]) is rewritten as 

y/^DtDpGixwrn) = ^(V* - ze^)(V M - ieAJG(x, y; m) 

= {dt-ieAt)y/^^{d v -ieA v )G{x,y-,m). (4) 

We want to expand Eq.(|l]) around R = 0. For this purpose we introduce the Riemann 
normal coordinates [ffHl . In the Riemann normal coordinates framework the metric tensor 
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is expanded as 



g(y) = -i- l -RfKv-xY{y-xy, 



(5) 



where the suffix (0) for R^L V and RffJ designates the curvature tensor at the origin x. 
Substituting Eq.(|) into Eq.(J|) we obtain 



^D^Dfiix^m) 



iT + -R {0) a p(y-x) a {y-xf^ 
o 

~R®\{y - x) a } (9„ - ieA v )G{x, y; m) 



(6) 



where we keep only terms independent of curvature or linear in curvature. We expand 
the GF, G(x,y;m), as 



G(x, y; m) = G (0) {x, y;m) + (x, y; m) + 0(R 2 



(7) 



where G^ and G^ represent the terms independent of R and the terms linear in R 
respectively. 

Substituting Eqs.(||) and (0) into ([!]) we can perturbatively solve the Klein-Gordon 
equation (|I]) about the spacetime curvature R. The piece independent of curvature in 
Eq.([l|) is given by 

- ieAjfa - ieA v ) + m 2 ] G {0 \x, y; m) = -5 D (x - y). (8) 

Thus G^^x, y \ m) satisfies the Klein-Gordon equation in flat spacetime. We will give the 
explicit expression for G^ \x, y; m) below. The piece linear in R in Eq.(|T]) becomes 



6 



R {a) a p{y-xY{y-xy^ 



--R^ a %{y-xY{y-xf 



(df, - ieA /Ji )(d v - ieA v ) 
-\R {0) \(y ' xY(du ~ ieA v ) + G®(x, y; m) 



+ 



rTidn - ieAJidv - ieA u ) + m 2 



x,y\m) 



0. 



(9) 



Hence the piece G^ involving the terms linear in R is expressed with the help of the GF 
G^(x,y;m) in flat spacetime (for more details of local momentum representation, see 

gHP|0). 
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It is more convenient to introduce the new variable defined by 

G(x, y; m) = y)G{x -y;m), 

where Q(x,y) satisfies 

{d^-ieA^{x,y) = 0. 
Inserting Eq.(|l(]) into Eq.(||) we obtain 

1 



(10) 
(11) 
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R {0) a p(y-xT{y-xfd^d, 



-3^ (0) a fi(v - x) a (v - x)%d> 
-lR {0)U a(y ~ x) a d u + f flH &°\y - x; m) 



+ (d,d, + m 2 ) G^ix-y 



-y:m) = 0. 



(12) 



^-dependence disappears in the relationship between &°^(y — x; m) and G^'(x — y; m). 
Therefore the GF & l \x — y; m) is given by 



& l \x — y; m) 



d D zGqq(x — z; m) 



- 2 -R {0)V a {z - y) a d u + ei? (0) ] &°\z - y- m) 

where the function G 00 (x — y;m) satisfies 

(d^ + m 2 )G 00 (x -y;m) = -5 D (x - y). 



(13) 



(14) 



The linear curvature correction terms of the scalar GF are expressed by the GF in flat 
spacetime. 

In the constant curvature spacetime, 



the equation ( |I~3"D simplifies to 

& l \x — y; m) 
R 



x 



(15) 



(16) 



D(D -I) J dDzG °°( x ~ Z] m ) 

D 3 Va(s(z - y) a {z - yfd^d, + \{z - y)»(z - yfd^ 
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--(D-l){z-yydp + D{D-l)Z 



G (0) {z-y;m). 



(17) 



As is shown in Eq.(§), G^°\x,y;m) satisfies the Klein-Gordon equation in flat space- 
time. According to the Schwinger proper-time method the Klein-Gordon equation is 
exactly solvable in the constant electromagnetic field [18, TTJ]| To calculate G^°\x, y; m) 
we introduce the proper-time Hamiltonian which is defined by 



H{x^id^\ 



[id 11 + eA»)(idp + eA^) - m 2 ] G (0) ( 



x,y;m) 



For this proper-time Hamiltonian the time evolution operator U(x,y; s) is defined by 



d 

i-^U(x,y; s) = H(x,id^)U(x,y; s), 



with the boundary conditions 

f lim s ^ U(x, y; s) 
\ lim s ^„ 00 U(x, y; s) 

Comparing Eq . (|19"D with Eq.(|TJ) we find 

G^°\x, y; m) = —i 



S D (x-y), 
0. 



dsll (x, y; s). 



(19) 



(20) 



(21) 



The time evolution operator U(x, y; s) is obtained by solving the equation of motion 
for x^ and 7r M = id^ — eA 1 . For a constant electromagnetic field, x M (s) = W (s)x fl U(s) 
and 7r M (s) = W (s)7r M t/(s) satisfy the following equations : 



dx»(s) 
ds 



dir^(s) 
ds 



i[H, x"} 



2eF^7r l 



These differential equations can be easily solved, 

7r( s ) = e 2eFs 7r(0), 



JleFs 



X(S) 



x{0) 



eF 



-vr(O), 



(22) 
(23) 

(24) 
(25) 



where we use the matrix notation. Substituting Eqs.(p4[) and ( |25| ) into Eq.(|T9"D we obtain 
the differential equation for the time evolution operator U(x,y; s), 



d 

i— U(x, y; s) = (vr 2 (s) - m 2 )U(x, y; s) 



|(x — y)K(x — y) — -tr[eFcoth(eFs)] — m 2 | U(x, y\ s), 



(26) 



^ For other types of exact solutions of the Klein-Gordon equation in the external electromagnetic 
field, see M. 
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where K is defined by 



K = ie 2 F 2 [sinh(eFs)]- 2 . 



(27) 



After the integration over s in Eq.(^) the time evolution operator is found to be : 



U(x,y; s) 



(47T 



,D/2 



Q(x,y)s D ^ 2 exp< — trln 



sinh(eFs) 
eFs 



x exp 



— (x — y)eF coth(eFs)(x — y) + im 2 s 



{21 



Substituting this time evolution operator (p8|) to Eq.(pT|) we find the GF, G^(x,y; 



m) 



1 

- - (47r) fl / 2 



dsU (x, y; s) 



/u 
dss~ D ' 2 exp 
-co 



1 , 

-~tr In 
2 



/ sinh(eFs)\ 
V eF~s ) 



x exp 



[x — y)eF coth(eFs) (x — y) + im 2 s 



Taking the limit F — » and changing the variable s —>■ — s in Eq.( 

G 00 (x -y;m), 



(29) 

we get the GF, 



G 00 (x -y;m) = - 
x exp 



c/.S' 



o (4tts) d / 2 



7T 2 

i-L> - —(x - y) M (x - y)^ - im 2 s 



(30) 



where the phase factor is determined to satisfy the boundary condition (|20|). 

In the present paper we consider the constant magnetic field along the z-axis. For this 
constant magnetic field, F 12 = — F 2 i = B, the GF, &°'(x — y;m), reduces to 



&°\x-y;m) ■■ 

eBs 
x ■ / p x exp 



ds 



o (4vrs) fl / 2 
-ijD - —(x - yyC^x - yf - im 2 s 



where 



= T)^ + F/F; 



. ; , 1 — eBs cot (eFs) 



B 2 



(31) 



(32) 



Inserting the Eqs.(|3!lD and ( PT| ) into ( ]TT| ) we obtain the GF linear in R 



& x '(x — y;m 
x exp 



R 



D(D-l) 



d D z 



dt 



ds 



4 4r ; v 7M 



o (47rt) D / 2 Jo (4^5)^/2 

2j 
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X 



D-3 
6 



Vapiz-yYiz-yfd^ 



+\{z - vY{z - vYd,d v - jj(-D - l)(z - y)"0„ + £>(D - 1)£ 



e5s 

X — — ; - ; exp 



sin(e5s) 



4 4s 



(33) 



Thus, we obtain the explicit expression of the scalar GF in the external gravitational and 
magnetic field. 

We need only the coincidence limit x — > y of the GF to calculate the effective potential 
of the models considered in the present paper. After the Wick rotation s — > —is,t — > 
— it, z° — > —iz° and the integration over z SpG(x,x;m) simplifies to 



SpG<°>( 



x,x;m) 



(4tt)^/ 2 Jo 



; -D 2 e ^ S ( 2 \ 

«SS ' ; ; - ; exp I — 771 Sj 



[x, x; m) 



R 



x 



(4tt) d / 2 D(D - 1) Jo Jo 

( s+t )(2-D)/2 g£ > 



sinh(e5s) 

(it (is exp — m 2 (s + t) 



(34) 



oo /-oo 



X 



+ 



1 + eBtcoth(eBs) sinh(el?s) 
D-3 
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D-3 



(£>-2 + 2e£scoth(e£s)) + 



-2D + 1 



(D - 2)t 

s + t 



6 s 



(D - 2 + 2e5s coth(e-Bs)) + 2j ° + - e£ coth(eBs) 



2t 



x 



1 + eBt coth(eSs) 
2(D - 1 



D{D-2) ( t \ 2 ( eBtcoth(eBs) \ 

' - eBt coth(eBs) J 



+ 



4 Vs + t/ ^1 + 

(D - 3) ((e5s) 2 coth 2 (e£s) + l) + AeBs coth(eBs 
t 



3s 

(D - 2)t 



s + t 1 + eBt coth(eBs] 



+ D(D-1)£ 



(35) 



Eqs.(|34D and fl35|) correspond to the the vacuum self-energy of the free scalar field with 
mass m at the one loop level. 



2.2 Green function for a spinor field 

Next we construct the spinor GF in an external electromagnetic and gravitational field. 
Let us write the GF, which obeys the Dirac equation: 

(ij^D^ - m)S(x,y;m) = — F =5 D (x-y), 



'-9 



(36) 



8 



where the covariant derivative D p includes the electromagnetic potential A 



i-i- 



1 

■M + 2 " 



(37) 



The local Dirac matrices 7 M (x) are expressed through the usual flat ones 7 a and tetrads 

(38) 



v 



Y(x) = 7 a e^(x), 

<?ab = T [7a, 7b]- 



The spin-connection has the form : 

1 
2 



= - d v e\) + - e ™e b »e cfl (d p et - d v e a p ) 

1 . _ 1 



-e bu {d,e a v - d v el) - - e b ^e c ^d p el - d v e c p ). 



(39) 



2 s-H>-v ~"~H/ 

Dimensions of the spinor representation are supposed to be four. Greek and Latin indices 
correspond to the curved and flat tangent spacetimes. 

To calculate the linear curvature corrections the local momentum expansion formalism 
as in the previous subsection is the most convenient one 15, T7j. In the Riemann 
normal coordinates framework the tetrads e At a (x) and the spin connection u ab fl a a b are 
expanded as 

1 6 (40) 
u%a ab = -R {0)a \x{y - x) x a ab . 

Substituting fl3"9|) and ([H]) into the (|3T|), we obtain the following equation for the GF : 



* 7 a (V a + \r^ p Uv - x) p (y - *y) (d„ + \R ( b % x (y - *) x ° bc - 

—m S(x,y;m) = S D (x — y). (41) 

Fulfilling the expansion on the spacetime curvature degrees 

S = S (0) + S w + • • • , (42) 

where is the GF in the flat spacetime, ~ 0(R) and so on, we receive the iterative 
sequence of equations: 



i + e -4-(x) — m 



S(°\x, y; m) = 5 D (x — y) 



(43) 



(i + e A"(x) — m)S ( - 1 \x, y; m) 

+ l -^R^ paa (y - X )P(y - xY{d, - ieA^x)) 
o 



S {0) {x,y;m) = 0. 



(44) 
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Here and below we can forget about the difference between the two kinds of indices 
(Greek and Latin) because it lies beyond of linear curvature approximation. 



We assume that just as in the flat spacetime GF has the form [18 

S(x,y;m) = $(x,y)S(x - y;m), 



(45) 



where the function $(x, y) is introduced in Eq. ([□]). Then, we can find the equation which 
determines the S^\x — y;m) function, excluding the evident dependence on A^{x): 

(i $ — m)S^ (x — y; m) 



- ---fR®* ( y - x)"{y - xYd^°\x - y- m) 
o 

\l a a hc Rf caX {y-xfS^{x-y-m). 



(46) 



The flat spacetime GF in the external electromagnetic field is supposed to be known 
TTj |l^ , 0. Denoting as S 00 (x — y; m) the GF, satisfying the equation 



(i ft - m)S o(x -y;m) = 5 D (x - y), 



we obtain 



J d D zS q(x — z; m)S ( - 1 \z — y;m) 

= -\i a R {0 KUx - y) p {* - y) a d^\x - y- m) 

-\l a a bc Rtl x {x-y) x ~SV\x- y] rn\ 



(47) 



(4* 



or, finally, 



(x — y;m) = J d D zS 00 (x — z;m) 
j a R^ a {z - y) f \z - yfd^iz - y- m) 
S (z - y;m). 



x 



% 
6 



-\r° bc RfL{z-y) x 



(49) 



However, we need only the coincidence limit x — > y to calculate the effective potential. 
It provides us the opportunity to simplify fl49"|) especially for the constant curvature 
spacetimes QT5p. Thus the expression for the GF, in the spacetime with an arbitrary 
dimension D is the following: 



S«(0; 



iR 



m) 



d D zSqq{—z; m) 



2^d^S {0 \z]m) 



12D(D - 1 

— 2z 2 "f "d„£< 0) (z; to) + 3(D - 1) ^S (0) (z; m) 



(50) 
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Now we begin the calculation of the Green function for the D = 4 Gross-Neveu model 
in the external magnetic field. According to the Schwinger proper-time method as in the 
previous subsection the flat spacetime GF is found to be [11, IS 



£(0) ( 



z: m = — i 



ds 



exp 



7T 

— i I — h sm' 
1 2 



exp 



As 



es 



Jo I6(irs) 2 

where C Mi/ is defined in Eq.(^). For Soo( — z\ m) we have directly from (|5ll 



(eSs) cot(eBs) - y^T"^ 



5 00 (-^; m) = -i 



exp 



. / 7T „ ■ 

-i - + m 2 t H 

\ 2 4t. 



m 

2t 



(51) 



(52) 



Substituting ( j52|) in 

SpS (1) (0;m) = 



16(7Tt) 2 

i and calculating the trace over the spinor indices, we have : 
iRm r d 4 zdtds 



36 

x exp 



[IQlxHs) 2 

. . 9 t + s 9 9 1 + e£>t cot(eSs) 
(i + s)m 2 + — z 2 - zl 



4ts 



If 



x<^ eSscot(eSs) 



2t 



9 f 9 A-7eBscot(eBs) 



z 2 ( 9 , 7 -4eBscot(eBs) 



2s 



2s 2 



zn 20.(1 - eSscot^s))' 



-e 2 5 2 



s ( 2z^ + -zj 



+-z 2 zj_(eBscot(eBs) - 1) 



where 



2,2 2 
^1 ~r ^2 5 Z \\ 



Z Q z 3 . 



(53) 
(54) 



After the Wick rotation and integration over z, one gets: 



SpS (1) (0;m) 



iRm 



dtds 



exp[— {t + s)m 2 



96vr 2 J (t + s) 2 (l + eBt coth(eBs)) 
x eBt(eBt + e£s) + 2(eBt + 3e£s) coth(eSs) 

+2eBt(eBs - eBt) coth 2 (e.Bs) 



(55) 



The same program can be done for the 3D case. The GF in the flat space-time with 
const; 

S i0) (z; m) 



the external constant magnetic field has the following form []TT| , |T2 

(is 



-i(7r/4+sm 2 ) 



(tts) 3 / 2 



exp(-— 
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es 



eBscotieBs)-—^^^ 



For S'oo(— z;m) we have directly from 

roc 

S 00 (-z;m) = -i 



8(vrt) 3 / 2 

Substituting in (0), we have: 



cxp 



m — 



2t 



SpS w {0;m) 



Rm r d 3 zdtds 



exp 



1152^7 (ts) 3 / 2 

2 1 + eBtcot(eBs) 
~ Zl it 

x <^ eBs cot(eBs) Zq(- 



3 3 - 2e J Bscot(e J Bs) 
t s 



2 .3 1 -2eBs cot( eBs) ^ 

-^zlziil - eBscot(eBs)) 2 
2s z 

+e 2 B 2 



-2s(zl + z 2 



+ -ZqZ 2 ± (1 - eBscot(eBs)) 
After the Wick rotation and integration over z, one gets: 



iRm 



dtds 



x 2eBt(eBt + eBs) + (9eBs + heBt) coth(eBs 
+eBt(eBs - 3eBt) coth 2 (eBs) 



exp[— [t + s)m 2 



(56) 



(57) 



(51 



(59) 



The proper-time integrations remain in our final expressions of the scalar and spinor 
Green functions in an external magnetic and gravitational field. All remained integrands 
are exponentially suppressed at the limit s — > oo and/or t — > oo. There are divergences 
at s -> and/or t -> for D > 2. 



3 NJL model 

Let us discuss now the NJL model in an external magnetic and gravitational field 
(for a review, see |J). The NJL model is one of the simplest models where dynamical 
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symmetry breaking is possible. It is well-know that the chiral symmetry of the model is 
dynamically broken down for a sufficiently large coupling constant. In the present paper 
we want to know the effect of an external electromagnetic and gravitational field to the 
dynamical symmetry breaking. 

The NJL model is defined by the action which is given by 

A 

where N is the number of the fermions. This action has the chiral U(l) symmetry. 
Introducing the auxiliary fields 



S 



{^) 2 + (ipij 5 ipf 



(60) 



we can rewrite the action (BOf) as: 



S 



d D x\ / —g 



N 



iipj^D^ip - ^(c 2 + vr 2 ) - ip(a + wr7 5 )V> 



(61) 



(62) 



In the leading order of the 1/N expansion the effective action is given by 
1 „ , , 



N' 



eff 



a, 7r 



d D x\J—g 



a 2 + ir 2 



2A 



i\ndet[i^(x)D^ - (a + i^n)]. (63) 



We can put now it = because the final expression would depend on the combination 
a 2 + 7r 2 only. In this case the effective action of the NJL model has the same form as the 
one of the Gross-Neveu model 0. Note that Eq.(|63"D represents the particular example 
of the effective action for composite fields [[HI] , 

Defining the effective potential as V e ff = 
tions of a and tt, we obtain 



T e ff/N J d XsJ—g for constant configura- 



a 



Veff(a) = — +tSp\n(x\[i^(x)D fl -a]\x). 



(64) 



The second term of the right hand side in Eq.flBip is rewritten as 

iSp\n(x\[i'y tJ "(x)D l j, — <r]\x) = —iSp / dm S(x,x;m). 



(65) 

Thus we can express the effective potential of the NJL model by the spinor Green function 
S(x, x, a). Substituting the Green function evaluated in the previous section to Eq.(|65|) 
we obtain the effective potential in an external magnetic and gravitational field 

In four-dimensional spacetime the linear curvature correction for the effective potential 
is given by 



'eff\ 



R 



dtds 



192vr 2 A/a* A/A" (t + s) 3 (l + eBtcoth(eBs)) 2 
x exp[-(t + s)a 2 }[eBt(eBt + eBs) + 2(eBt + 3eBs) coth(eBs) 
+2eBt(eBs - eBt) coth 2 (eBs)}, (66) 
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Figure 1: The behavior of and V^A are shown with the varying magnetic field 
eB(= 0, /i 2 /2, /i 2 , 3/i 2 /2, 2/i 2 ) for fixed A(= 1/2.5) and fixed A(= 10/x) in four dimensional 
flat spacetime. \x is an arbitrary mass scale. We normalize that V e f f (0) = V e f f (0) = 0. 



where we introduce the proper-time cut-off A. In four dimensions four-fermion models 
are not renormalizable. Here we define the finite theory by introducing the proper-time 
cut-off A. Taking into account the proper-time representation for V e ff in four- dimensional 
flat spacetime ]n|, we obtain the effective potential with the linear curvature accuracy: 

cr 2 1 f°° ds 
Veff{°) = ^x + ^2 J i/A2 ^exp(-sa 2 )eBscoth(eBs) 

R r 00 r°° dsdt 



192vr 2 Ji/A* Ji/a» (t + s) 3 (l + eBt coth(eBs)) 2 
x exp[-(t + s)a 2 }[eBt(eBt + eBs) + 2(eBt + 3eBs) coth(eBs) 
+2eBt{eBs - eBt) coth 2 (eBs)). (67) 

It should be emphasised here that for B = we obtain just the same expression for 4D 
NJL model V e ff in the proper-time representation as it has been found out in |TJ| (See 
also §). 

For D = 3 the linear curvature correction for the effective potential is given by 

r( i) , R r°° r°° dtds 



eff ^ ~ ~ 144vr 3 / 2 Ji/a» Ji/a? (t + s) 5 / 2 (l + eBtcoth(eBs)) 2 
x exp[-(t + s)a 2 }[2eBt(eBt + eBs) + (9eBs + 5eBt) coth(eBs) 
+eBt{eBs - 3eBt) coth 2 {eBs)]. (68) 
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Veff(cT)/fi 4 




-0.008 - 



0.2 0.4 0.6 OA 

a/fi 

(a) V eff for R = 3/j 2 . 



-0.002 - 




-0.008 - 







0.2 0.4 0.6 0.1 

a/fi 

(b) V eff for eB = fi 2 /2. 



Figure 2: The behavior of the effective potential V e ff is shown with the varying B or R 
for fixed A(= 1/2.5) and fixed A(= 10//) in four dimensions. 



For D = 3, four-fermion models are known to be renormalizable in the sense of the 1/N 
expansion. Thus the cut-off dependence in the effective potential disappears after the 
usual renormalization procedure. Taking into account the proper-time representation 
for V e ff in flat spacetime (Tl|, |12|| , we can write the effective potential with the linear 
curvature accuracy: 



Veff (a) 
R 



a 



+ 



2A 4vr 3 /2 J 1/A2 s 5/2 



ds 



exp(—sa 2 )eBs coth(eBs) 
dsdt 



144^3/2 J 1/A2 J 1/A2 ( t + s )5/2(! + e Btcoth{eBs)) 2 



x exp[— (t + s)cr 2 



2eBt(eBt + eBs) + (9eBs + 5eBt) coth(eBs) 



+eBt(eBs - 3eBt) coth 2 (eBs) 



(69) 



The effective potential ([ 
in the proper-time representation with B 



gives again the correct expression for V e ff of 3D NJL model 
JR. 

Using Eq. ( |6?| ) we numerically calculate the effective potential for D = 4 and show 



it in Figs.l and 2. In drawing Fig.l and 2 the coupling constant A is kept in the region 

where the chiral symmetry is broken down for R = B = 0. To see the curvature effect 

in the external constant magnetic field we divide the effective potential into the terms 

-(1) 



independent of the curvature and the terms linear in the curvature, K// — V e ff + K 



In Fig.l typical behaviors of the effective potential Vgfl and V^i are given for fixed four- 



eff 



eff 



fermion coupling constant A and proper-time cut-off A. It is clearly seen in Figs.l and 
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2. (a) that a constant magnetic field B decreases the potential energy in the true vacuum 
even in the curved spacetime. Curvature effects depend on the sign of the spacetime 
curvature R as is shown in Fig. 2.(b). A negative curvature decreases the potential 
energy. On the contrary, positive curvature increases the potential energy. The broken 
chiral symmetry is restored for a sufficiently large and positive curvature even in an 
external large magnetic field. 



4 SUSY NJL model 



In the present section we consider the supersymmetric NJL model in the external constant 
magnetic field and external gravity. 

In flat space such a model has been introduced in Ref. \12\. Generalization of the 
model for the case of non-minimal coupling with the external gravitational field (via 



scalar-gravitational coupling constants) has been presented in Ref. p3| . We consider the 



model of Ref.|23[ in curved spacetime with the external magnetic field. The action in 
components is given by 



S 



-0 f + a 2 + ^R)<f> - ct (ZW^ + a 2 + £ 2 i?)0 c 

1 



2A 



a 



(70) 



where a is an auxiliary scalar as in the original NJL model, ip is N component Dirac 
spinor. 

Note that actually the action (|70| ) represents SUSY NJL model non-minimally inter- 
acting with the external gravity and minimally interacting with the external magnetic 
field. It is evident that supersymmetry of SUSY NJL model is always broken in the 
external fields. So it maybe natural to call the theory with action ( |T0"D the extended NJL 
model. 

At the leading order of the 1/N expansion the effective potential of the SUSY NJL 
model (|70|) is given by 

l -a 2 + iSp\n(x\[i^{x)D tl - a)\x) 



eff 



(J 



2A 



-iSplnfcl [£>"£>„ + a 2 + ^R}\x) 

-iSp\n(x\[D»D^ + o 2 + &R]\x). (71) 

The second term in the right hand side of Eq. (|7T|) corresponds to the radiative correction 
by spinor fields and is rewritten by the spinor GF as is evaluated in the previous section 
(See Eq.(^)). The third and the fourth terms in the right hand side of Eq.flTTD correspond 
to the radiative correction by scalar fields and is represented by the scalar GF 



iSp ln(x| [D^D^ + a 2 + £R] \x) 
2i I m dm G(x, x, m). 



(72) 
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Because of the supersymmetry both the radiative correction by spinor fields and the one 
by scalar fields are cancelled out in flat spacetime with vanishing electromagnetic fields. 
Thus the vacuum expectation value of a disappears for R = B = 0. Substituting Eqs.( 
and (|35D into Eq.(^) and integrating over m one gets 

2i m dm G(x,x,m) 



(4tt 
1 



1 f°° j -DI2 eB ( 2 \ 

Wo dss ^Bs) e ^y- as ) 



R 



oo roo 



(4tt) d / 2 D(D-1)Jo Jo 

{ s + ty D ' 2 



sinh (eBs 

dt ds exp — (s + t)a 2 

eB 



1 + eBtcoih(eBs) sinh(e5s) 
D - 3 



6 

D - 3 



(D-2 + 2eBs coth(eBs)) + 



-2D + 1 



(D - 2)t 

s + t 



6 s 



(D-2 + 2eBs coth(eSs)) + 2D + - eB coth(eBs) 



2t 



x ■ 



1 + eBt coth(eBs) 



2(D-1] 



D(D - 2) ( t 



4 



s + t 



+ 2 



eBtcoth(eBs) 
1 + eBt coth(eBs) 



1 



x 



(D - 3) ((eBs) 2 coth 2 (eBs) + l) + 4eEs coth(e5s; 



3s 

(£> - 2)t 



s + t 1 + eBtcoth(eBs) 



D{D-l)i 



(73) 



In some special dimensions Eq.([73|) simplifies. Inserting the Eq.([T3|) into Eq.([7T|) we 
obtain the final expression of the effective potential of the SUSY NJL model. Therefore 
the effective potential for D = 4 reads 



Veff(a) 
R 



oo roo 



dsdt 



192vr 2 A/A" A/A2 (s + t) 3 (l + eBt coth(eBs)) 2 
x exp[-(s + t)a 2 ][eBt(eBt + eBs) 

+2(eBt + 3eBs) coth(eBs) + 2eBt(eBs - eBt) coth 2 (eSs)] 



00 ds eBs 



8n 2 A/a 2 s 3 sinh(ei?s) 
R 



exp 



oo roo 



dsdt 



48?r 2 A/a 2 A/a 2 (s + t) 2 (l + eBt coth(eBs)) 
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X 



eB 



+ 



sinh(e-Bs) 
1 4 



6 s 



exp[— (s + t)o A 
eB coth(e.Bs) 



4 1 

eBs coth(e£>s) 

3 6 

2t 



2t 
s + t 



eBt coth(eBs) 



+2 



t 



eBtcoth(eBs) 



+ 



3s 



t 



\l + eBt coth(eBs) 

(eBs) 2 coth 2 (eBs) + 1 + AeBs coth(eBs) 
t 



s + t 1 + eBt coth(eBs) 



+ 3(6 + 6 



(74) 



where we regularize the effective potential by introducing the proper-time cut-off A. The 
effective potential (]74D reproduces the expression for V e ff(a) in Ref.p3| at the limit 
B — > 0. The proper-time integrations in Eq.([74"D are numerically performed below. 
For D = 3 the effective potential reads 



VeffW 

R 



ds 



2A 4vr 3 / 2 A/a? s 5 / 2 



a 

+ 



exp(— sa 2 )eBs coth(eSs) 
dsdt 



144tt 3 / 2 A/A2 Ji/A* (t + s) 5 / 2 (l + eBtcoth(eBs)) 2 



x exp[— (t + s)a 2 



2eBt(eBt + eBs) 



-(9eBs + 5eBt) coth(eBs) + eBt(eBs - 3eBt) coth 2 (eBs) 



1 



00 ds 



eBs 



4tt 3 / 2 Ji/a' s 5 / 2 
R 



exp( — sct 



smh(eBs) 
dsdt 



24tt 3 / 2 Ji/a' h/tf (t + s) 3 / 2 (l + eBtcoth(eBs)) 



x 



+ 



eB 



sinh(ei?s 
t ^ 2 



exp[— (s + t)a z 



5 t 5 2eBt coth(eBs) 



3 s + t 3 1 + eBt coth(eBs) 



+ 



4 t 



8 / e5tcoth(e J Bs) \ 2 
3 V 1 + eBt coth(eBs) ) 
eBtcoth(eBs) 



3 s + 1 1 + eBt coth(eBs) 



(75) 



Since the contribution of scalar fields and fermion fields are canceled in a supersymmet- 
ric theory only a symmetric phase is realized for R = B = 0. The external gravitational 
and magnetic fields considered here break the supersymmetry of the theory. Thus we 
expect that the chiral symmetry may be broken down. Using Eq. (|74|) we numerically 
calculate the effective potential of the SUSY NJL model for D = 4 to study the phase 
structure of the theory. 
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K ( /V)/^ 
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(a) V}jj for R = 0. 



0.012 








eB = 



0.2 0.4 0.6 0.8 



(b) for = and 6 + 6 = 0. 




-0.008 



0.2 0.4 0.6 0.8 

(c) V e f f for R = and & + £ 2 = 1. 



Figure 3: The behavior of the effective potential V e jj- and V^j are shown with the varying 
magnetic field eB{= 0,5/x 2 , 10/i 2 , 15/x 2 ) for fixed A(= 1/2.5), and fixed A(= 10//) in four 
dimensions. // is an arbitrary mass scale. We normalize that V e ff(0) = 0. 
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(a) V eff for R = 5/i 2 and 6 + 6 = 0. 
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(b) V eff for R = 5/i 2 and & + £> = 1. 
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0.1 

(c) for e_B = 15/x 2 and £i + £ 2 



0.2 
= 0. 




(d) V eff for eB = 15/i 2 and & + & = 1. 



Figure 4: The behavior of the effective potential V e ff are shown with the varying B or 
i? for fixed A(= 1/2.5), and fixed A(= 10/x) in four dimensions. 
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In Figs. 3 and 4 we illustrate the typical behaviors of the effective potential for D = 4 
in the SUSY NJL model. In Fig.3 we draw the terms independent of R, V^°\ and the 
linear curvature correction to the effective potential, . The linear curvature correction 
strongly depends on the coupling constant £i + £2- There are large cancellations between 
the corrections of fermion and scalar fields. Thus the terms of the linear curvature 
correction in the SUSY NJL model are significally smaller than that in the NJL model. 
As is clearly seen in Fig. 4 the chiral symmetry is broken down for a sufficiently large 
magnetic field and/or a negative curvature for R = 5fi 2 or B = 15yU 2 . The external 
curvature has the opposite effects in the case of large £1 +^(> 1/2) for a small B and/or 
large a (See Fig.3. (c)). 

5 Conclusion 

In the present paper we discussed the combined effect of the gravitational and magnetic 
fields to the chiral symmetry breaking in NJL and SUSY NJL models. Chiral symmetry 
is broken at non- zero and non- positive curvature. On the same time, positive curvature 
acts against chiral symmetry breaking. Nevertheless, the magnetic field effects may be 
significally stronger in realistic situations corresponding to early Universe with primordial 
magnetic fields. 

It should be noted that the effective potential calculated in the present paper depends 
on how to introduce cut-off parameter in the proper-time integral, though we do not 
develop it any further here. 

Using the results presented in Ref. || it is not difficult to calculate the effective 
potential in above two models exactly on the constant curvature spacetimes (DeSitter 
or anti- DeSitter, for example || [7|). Such calculation shows that curvature effects 
estimated in this work are taken into account qualitatively correctly (one can consider 
then very strong curvature as well). However it is not clear how to make such calculation 
exactly both on curvature and on magnetic field. The only possibility to do so is to work 
on some gravitational- magnetic background where exact solutions of field equations are 
known (for example, on conformal spacetime with magnetic field). 

Another interesting proposal could be to study the dynamical symmetry breaking 
on electromagnetic- gravitational background representing a combination of constant 
electromagnetic field with gravitational wave (it is known the effective Lagrangian in 
pure electromagnetic or gravitational wave is trivial and no effect is expected). Such 
background may be roughly considered as signal which comes from strongly- graviting 
object with strong magnetic field. Then one may speculate on the possible use of the 
dynamical symmetry breaking in gravitational waves detectors. 
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